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An  anti-windup  controller  modification  is  implemented  in  control  system  design 
for  a  model  of  the  longitudinal  dynamics  of  an  air-breathing  hypersonic  vehicle. 
Anti-windup  control  allows  the  input  constraints  to  be  considered  explicitly  in  the 
design  of  linear  controllers  to  track  a  reference  trajectory  for  the  vehicle  velocity, 
altitude,  and  angle  of  attack.  The  presence  of  anti-windup  alleviates  the  need  of 
keeping  large  penalties  on  the  magnitude  of  the  control  input  to  avoid  the  occur¬ 
rence  of  saturation.  This,  in  turn,  allows  tighter  tuning  of  the  controller  gains  to 
obtain  faster  and  more  accurate  trajectory  tracking.  The  paper  employs  recent 
developments  in  anti-windup  design  to  deal  with  the  presence  of  exponentially  un¬ 
stable  dynamics,  which  are  typically  encountered  in  air-breathing  vehicle  models. 
Simulation  results  on  a  fully  nonlinear  model  are  presented  to  validate  the  con¬ 
troller  design. 

Keywords:  anti-windup,  tracking,  hypersonic,  air-breathing,  control 

I.  Introduction 

Seen  as  a  possible  way  of  making  space  access  more  affordable,  air-breathing  hypersonic  vehicles 
(HSVs)  have  been  studied  sporadically  over  the  last  five  decades  as  the  technical  advances  necessary 
for  their  development  were  missing.  New  technology  has  brought  a  renewed  interest  in  the  research 
area,  as  demonstrated  by  the  successful  flights  of  NASA’s  X-43A  in  April  and  November  of  2004. 
Air-breathing  hypersonic  vehicles  are  characterized  by  their  unique  design,  incorporating  a  super¬ 
sonic  combustion  ramjet  engine  located  beneath  the  fuselage.  This  esoteric  configuration  results  in 
strong  coupling  between  the  thrust  and  pitch  dynamics  of  the  vehicle,  which  in  combination  with 
flexible  effects  and  static  instability  make  the  vehicle  a  challenging  application  for  control.  Control 
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of  air-breathing  hypersonic  vehicles  is  a  difficult  task  because  of  the  strong  interaction  between 
the  airframe  and  the  engine  dynamics  resulting  in  sizable  couplings  between  thrust  and  pitch  mo¬ 
ment.^’ Furthermore,  due  to  the  vehicle  length,  relative  slenderness,  and  light  weight,  the  vehicle 
dynamics  is  characterized  by  the  presence  of  flexible  modes  which  interact  with  the  aerodynamic 
forces  and  moments. 

In  the  literature,  models  of  the  longitudinal  dynamics  of  air-breathing  hypersonic  vehicles  have 
been  developed,^’ with  a  few  concerted  efforts  to  incorporate  guidance  and  control.^’ 

In  the  cited  references,  mostly  linearized  models  of  the  vehicle  dynamics  have  been  considered  for 
control  design,  with  the  noticeable  exception  of  Tournes  et  al.,^^  where  nonlinear  variable  structure 
control  has  been  adopted.  On  the  other  hand,  for  hypersonic  vehicle  models  with  conventional 
actuation  and  negligible  flexibility,  the  use  of  nonlinear  control  system  design  methodologies  is 
more  common.^’ Like  many  real  world  systems,  air-breathing  hypersonic  vehicles  have  input 
saturations  which  must  be  considered  in  control  design.  Hard  constraints  on  the  inputs  can  lead  to 
performance  degradation  and  even  instability.  Anti-windup  control  provides  a  method  for  resisting 
this  performance  degradation  and  helping  to  avoid  the  onset  of  instability.  Several  anti-windup 
control  methods  are  discussed  in.^’^A6,i7 

In  our  previous  work,^  a  linear  set-point  tracking  controller  was  developed  for  a  particular  model 
of  the  longitudinal  dynamics  of  an  air-breathing  HSV  as  developed  in  Bolender  and  Doman.^  In 
that  work  however,  the  constraints  on  the  inputs  were  not  directly  taken  into  account.  Instead,  they 
were  accounted  for  manually  by  tuning  the  gains  to  keep  the  inputs  in  their  valid  operating  ranges. 
With  an  anti-windup  controller  the  constraints  on  the  inputs  can  be  taken  into  consideration  at 
the  level  of  the  controller  design.  In  addition,  by  re-tuning  the  gains  of  the  nominal  controller  and 
increasing  the  speed  of  the  reference  trajectories,  we  can  elicit  a  faster  response  from  the  system. 
Anti-windup  control  will  also  provide  some  measure  of  robustness  to  the  physical  constraints  on 
the  inputs  in  terms  of  stability. 

Designing  an  anti-windup  controller  for  an  air-breathing  hypersonic  vehicle  is  a  difficult  and 
challenging  task,  largely  because  the  linearized  model  of  the  air-breathing  HSV  is  unstable.  Most 
anti- windup  control  methods  only  deal  with  stable  plants  with  input  saturation  such  as  the  methods 
discussed  in.^  Global  stabilization  of  a  system  with  input  saturation  is  not  possible  for  systems 
that  are  exponentially  unstable  even  with  anti-windup  control.  In  the  latter  case,  the  role  of 
anti-windup  control  is  to  enlarge  the  region  of  attraction  for  the  equilibrium  of  the  system  with 
input  saturation.  The  work  by  Teel^^  provides  a  method  for  dealing  with  exponentially  unstable 
plants.  In  this  paper,  Teel’s  method  is  applied  to  a  case  study  on  linear  controller  design  for 
a  recently  developed  nonlinear  model  of  the  longitudinal  dynamics  of  an  air-breathing  HSV.  By 
means  of  computer  simulations,  it  is  shown  that  with  a  judicious  selection  of  the  gain  parameters, 
Teel’s  anti-windup  design  is  capable  of  enlarging  the  domain  of  stability  achieved  by  baseline  linear 
controllers,  and  significantly  improving  the  performance  of  the  closed-loop  system  in  terms  of  speed 
and  fidelity  of  response  in  tracking  setpoint  references  trajectories. 

The  paper  is  organized  as  follows:  In  Section  H,  we  provide  some  background  on  anti- windup 
design,  and  some  details  on  the  method  proposed  in^^  to  deal  with  unstable  plants.  In  Section  HI, 
we  briefly  introduce  the  model  of  the  longitudinal  dynamics  of  the  air-breathing  HSV  under  con¬ 
sideration,  and  preset  the  baseline  linear  controller  employed  in  the  case  study.  Then,  we  present 
the  implementation  and  discuss  the  results  achieved  by  the  anti-windup  modification  in  Section  IV. 
Finally,  we  draw  some  conclusions  in  Section  V. 

II.  Background 

Controller  windup  occurs  when  the  control  actuators  reach  their  physical  limits  or  saturations. 
This  results  in  a  difference  between  the  output  of  the  nominal  controller  and  the  actual  input  to  the 
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Figure  1.  Block  Diagram  of  Classical  Anti-windup  Modification 


plant.  Because  of  this  difference  the  states  of  the  system  are  incorrectly  updated  causing  controller 
windup.  Since  the  linear  controller  does  not  account  for  these  input  nonlinearities  performance 
degradation,  overshoot,  and  even  instability  can  result.  Many  ad  hoc  methods  have  been  developed 
over  the  years  to  modify  pre-designed  linear  controllers  to  account  for  input  saturation  without  local 
performance  degradation.  As  mentioned  previously,  several  of  these  methods  are  discussed  in.^  The 
majority  of  these  methods  require  the  plant  to  be  stable,  and  are  referred  to  as  ‘observer-based’  or 
‘classical’  anti-windup  scheme.  Consider  a  plant  model  given  in  the  form 

Xp  —  ApXp  -j—  BpS^ij(^Up^  (l) 

Up  —  CpXp 

with  state  Xp  G  control  input  Up  G  M^,  and  regulated  output  G  M^,  endowed  with  a  nominal 
controller  given  by  the  system 

Xq  —  AqXq  -\-  BqCp 
yc  =  CcXc  +  DcCp 

where  Xc  G  G  M^.  The  signal  Cp  G  denotes  the  tracking  error  Cp  —  yr  —  yp^  being  yr  a 

reference  trajectory.  A  vector- valued  saturation  nonlinearity  of  the  form 

^max  ii  U  ^  ^max 
U  if  Amin  ^  U  Amax 

Amin  if  U  ^  Amin 

is  present  at  the  plant  input,  where  the  inequalities  above  are  intended  to  hold  componentwise, 
and  Amax  and  Amin  are  vectors  of  upper  and  lower  limits,  respectively.  We  assume  that  under 
the  interconnection  Up  —  yc  the  closed  loop  system  without  input  constraints  is  well  posed  and 
internally  stable. 

The  ‘classical’  anti-windup  control  method  basically  reduces  to  an  output  injection  on  the 
nominal  controller.  The  nominal  controller  of  the  system  is  modified  to  be 

Xc  =  AcXc  + BcC  +  L[sdX{u)  —  u] 
u  —  CcXc  +  DcC  . 

Figure  1  shows  a  block  diagram  of  the  observer-based  scheme.  Methods  of  this  type  are  limited  in 
their  ability  to  control  unstable  systems:  in  deriving  the  error  system  between  the  constrained  and 
unconstrained  closed  loop  systems,  it  is  readily  seen  that  the  eigenvalues  of  the  plant  matrix  Ap  are 
part  of  the  set  of  eigenvalues  for  the  error  system.  Thus  if  Ap  has  unstable  eigenvalues,  then  the 
constrained  closed  loop  system  will  diverge  exponentially  from  its  unconstrained  counterpart  when 
the  inputs  are  saturated.  The  work  by  Teel  and  Kapoor^^  provides  the  basis  for  an  anti-windup 
control  method  for  unstable  systems.  In  the  cited  reference,  the  authors  provide  the  following 
definition  of  the  anti-windup  problem: 
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Figure  2.  Block  Diagram  of  Separate  Anti-windup  Controller  Method 


1.  To  have  the  closed  loop  system  be  unmodified  in  the  absence  of  input  saturation 

2.  To  provide  C2  stability  for  the  error  system  created  by  the  difference  between  the  system  with 
input  constraints  and  the  system  without  input  constraints 

The  controller  proposed  in^^  modifies  the  closed  loop  system  both  at  the  input  and  the  output  to 
the  nominal  controller  (see  Figure  2).  This  anti- windup  controller  is  given  by 

i  =  Api  +  Bp[s&X{yc  +  ^^l)  -  Vc] 

Vi  =  k{C}  (2) 

V2  =  —Cp^ 

where  is  a  feedback  term  that  stabilizes  the  error  system  (the  mismatch  between  the  constrained 
and  unconstrained  closed  loop  systems)  in  £2-  The  anti- windup  controller  changes  the  closed  loop 
system  by  modifying  the  interconnections  to  Up  —  sat(2/c  +  vi)^Uc  =  Vp  —  Vr  +  "^2-  Though  this  is 
still  an  observer  based  solution,  it  creates  n  additional  states  giving  enough  degrees  of  freedom  to 
blend  the  predesigned  linear  controller  with  a  nonlinear  anti- windup  controller. 


A.  Teel’s  Method  for  Anti-Windup  Design  for  Unstable  Linear  Systems 

The  problem  of  anti-windup  for  exponentially  unstable  linear  systems  is  far  more  difficult  than 
for  stable  systems.  In  unstable  systems,  once  all  inputs  are  saturated,  the  trajectories  may  grow 
unbounded.  The  goal  of  anti-windup  control  for  unstable  systems  is  to  enlarge  the  domain  of 
attraction  of  the  equilibrium.  The  paper  by  Teel^^  illustrates  these  difficulties  and  derives  a  solution. 
This  method  takes  the  anti-windup  solution  given  by  Teel  and  Kapoor^^  and  modifies  it  to  account 
for  the  unstable  modes.  To  this  end,  the  system  must  first  be  transformed  by  separating  the 
unstable  modes  from  the  stable  ones  such  that  the  system  can  be  given  in  the  form 


y 


f  Ai2 
\  0  A+ 


(3) 


where  and  represent  the  states  associated  with  the  stable  and  unstable  modes  of  the  system 
respectively.  Letting  x  =  {x^  fhe  system  (3)  is  written  as 


X  =  Ax  +  Bsa,t{u) 
y  =  Cx 


(4) 
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Figure  3.  Region  of  Attraction  for  the  Unstable  Mode  and  Operating  Region  for  /3(x+) 


where  A  =  T  ^ApT^  B  —  T  and  C  —  CpT.  The  solution  of  the  anti- windup  problem  for 

exponentially  unstable  systems  proposed  by  Teel  employs  a  controller  of  the  form 

i  ^  A^  +  B[sat{yc  +  ^;i)  -  yd 

vi  =  [P-l]yc  +  a{x+-P[x+-^+],PK{^))  (5) 

V2  =  -C^  -  D[sat{yc  +  vi) -yd 

where  /3(x+)  is  a  function  of  x^{t)  that  takes  values  between  zero  and  one,  and  a^(^)  is  a  feedback 
that  provides  C2  stability  of  the  error  system  governing  the  mismatch  between  the  constrained 
and  unconstrained  closed  loop  systems.  The  proof  is  given  in,^^  and  will  obviously  be  omitted 
here;  however,  to  understand  the  controller  behavior,  it  is  necessary  to  introduce  some  additional 
notation.  Let  T+  be  defined  as  the  region  of  the  state  space  in  which  the  state  corresponding 

to  the  unstable  mode,  can  be  driven  back  to  the  origin,  and  let  the  set  be  conservatively  chosen 
to  be  contained  in  T+  (see  Figure  3).  The  anti- windup  controller  given  by  equation  (5)  is  divided 
into  two  separate  controllers  that  span  three  different  operating  regions  defined  by  the  function 
/3(x+)  as  shown  in  Figure  3.  The  function  /3(x+)  is  a  design  parameter,  whose  choice  affects  the 
performance  of  the  anti-windup  controller.  When  G  then  /3(x+)  =  1.  When  leaves 
then  /3{x-^)  linearly  tapers  off  until  it  reaches  zero.  In  general  the  function  /3(x+)  can  be  chosen  to 
be  a  piecewise  linear  function  interpolated  and  extrapolated  from  chosen  data  points.  These  data 
points  are  related  to  the  regions  and  T+.  In  the  first  operating  region  {x^{t)  G  /3(x+)  =  1), 
the  anti-windup  controller  takes  a  form  similar  to  the  one  given  by  equation  (2)  written  as 

^  +  B[sat{yc  +  vi)  -  yd 

vi  =  a(^+,K(e) 

V2  =  -C^- D[sat{yc  +  vi) -yd- 

This  will  be  referred  to  as  the  nominal  anti-windup  controller,  and  it  can  be  designed  in  a  similar 
manner  as  the  anti-windup  controller  for  stable  systems  presented  in.^^  The  function 
can  be  chosen  as 

a(^+,«(e)  =  fcie^++«(e- 

The  feedback  gain  ki^  is  designed  to  make  Aj^  +  Bj^ki^  Hurwitz  using  LQR.  Weighting  is  specifically 
placed  on  the  unstable  state  to  keep  the  trajectory  bounded.  We  can  then  take  the  function 
to  be  the  linear  function  =  —B^P^  =  ~^2eC*  The  feedback  gain  ^2^  is  computed  optimally 
using  LQR  methodology.  This  guarantees  C2  stability  of  the  error  system  as  long  as  x^{t)  G 

The  system  response  may  require  that  x^{t)  leaves  the  region  In  this  case,  the  controller 
enters  the  second  operating  region,  and  the  anti-windup  controller  takes  the  form  given  by  equa¬ 
tion  (5).  In  this  transitionary  regime,  both  the  nominal  anti- windup  controller  and  an  auxiliary 
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anti-windup  controller  are  active,  with  /3(x+)  determining  the  contribution  from  each  controller. 
The  auxiliary  anti- windup  control  serves  the  purpose  of  keeping  x^{t)  G  T+  so  that  C2  stability  is 
maintained  and  has  the  form 


C  =  -4^  +  B[s&t{yc  +  fi)  -  Vc] 
vi  =  -yc  +  a{x+,0) 

V2  =  -C^- D[sat{yc  +  vi) -yc]. 

In  the  third  operating  region,  the  auxiliary  anti-windup  controller  takes  over  and  the  nominal  anti¬ 
windup  controller  is  turned  off.  If  goes  outside  of  T+  the  trajectories  will  grow  without  bound 
since  the  unstable  state  will  not  be  able  to  be  driven  back  to  the  origin  using  the  constrained 
input.  Thus  the  goal  of  the  auxiliary  anti- windup  controller  is  to  keep  G  T+  and  drive  the  state 
back  into  the  region  Obviously,  anytime  the  auxiliary  anti- windup  controller  is  used,  the 
response  of  the  system  will  be  altered  more  than  for  the  nominal  anti-windup  controller.  This  is 
because  the  goal  of  the  overall  control  scheme  has  shifted  from  trying  to  track  the  reference  input, 
to  trying  to  keep  the  unstable  state  from  diverging.  In  this  case,  the  reference  trajectory  which  is 
commanded  is  not  achievable  using  the  constrained  control  inputs,  so  the  system  response  deviates 
from  the  reference  trajectory. 

III.  Anti- windup  Implementation  for  Air-breathing  HSVs 

The  longitudinal  dynamics  for  the  air-breathing  hypersonic  vehicle  model  considered  in  this 
paper,  as  given  in  Bolender  and  Doman,^  is  described  by  the  following  nonlinear  equations 


Vt  =  —  cos  a  —  D)  —  g  sm{6  —  a) 


a  = 

(  Tsina  L)  +  Q  +  cos 
mVt  Vt 

(9  —  a) 

lyyQ  — 

M  +  Xpifii  +  ^p2V2 

h  = 

Vt  sin(0  —  a) 

9  = 

Q 

kiiii  = 

-2CioJirii  Ni 

1 

1 

^2^im 

dyy 

k2m  = 

—  2C,2^2'n2  —  ^292  +  N2 

1 

1 

^1^292 

^yy 

where  T^L^D^Ni  are  the  thrust,  lift,  drag,  and  generalized  elastic  forces  respectively,  M  is  the 
pitching  moment  about  the  ^-axis,  lyy  is  the  moment  of  inertia,  uJi  are  the  damping  coefficients 
and  natural  frequencies  of  the  elastic  modes.  Furthermore, 

ki  =  l  +  p- ,  i  =  1, 2  =  /  rhf^(j)f{^)d^  ^2=  [  ma^(pa{Od^  > 

^yy  Jo 

where  rhf^  iha  are  the  mass  densities  of  the  forebody  and  aftbody  respectively  and  (pa  are  the 
mode  shapes  for  the  forebody  and  aftbody  respectively.  The  reader  is  referred  to  Bolender  and 
Doman^  for  further  details. 

A  linearized  model  about  a  trim  condition  of  the  nonlinear  vehicle  dynamics  described  by  the 
equations  above  is  given  in  the  form  (1),  where  Xp  G  M^,  Up  G  M^,  and  yp  G  represent  deviations 
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States 

Control  Inputs 

Vt 

Vehicle  velocity 

4 

Control  surface  deflection 

a 

Angle-of-attack 

ATo 

Total  temperature  change  across  combustor 

Q 

Pitch  rate 

Ad 

Diffuser  area  ratio 

h 

Altitude 

Controlled  Outputs 

e 

Pitch  angle 

Vt 

Vehicle  velocity 

Vi 

Generalized  elastic  coordinates 

a 

Angle-of-attack 

Vi 

Time  derivative  of  the  elastic  coordinates 

h 

Altitude 

Table  1.  Definitions  of  the  states  and  inputs 


States 

Inputs 

[ft/s] 

7846.4 

6  [rad] 

0.13716 

a  [rad] 

0.034217 

ATo  [deg  R] 

484.49 

Q  [rad/s] 

0 

0.35005 

h  [ft] 

85000 

6  [rad] 

0.034217 

m 

1.6105 

m 

1.4582 

Vi 

0 

Table  2.  Trim  condition 


of  the  corresponding  trajectories  of  the  non-linear  system  from  the  trim  condition.  States,  inputs 
and  outputs  are  arranged  in  the  same  order  as  they  are  given  in  Table  1,  while  the  trim  condition 
is  given  in  Table  2.  The  final  performance  objective  is  to  achieve  a  simultaneous  setpoint  tracking 
of  the  angle  of  attack  ((a),  velocity  (1^),  and  altitude  (h)  without  the  actuators  violating  specific 
constraints,  given  in  Table  3.  The  plant  has  been  augmented  with  integrators  to  remove  the 
presence  of  steady-state  error  in  setpoint  tracking.  After  integral  augmentation,  the  model  is 
written  as  in  (4).  The  nominal  controller,  based  on  the  synthesis  given  in,^  is  a  linear  optimal 
tracking  controller  that  minimizes  the  cost  function 

J  =  -  /  (e^Qe  +  u^Ru)  dt^ 

2  Jo 

where  Q  =  >  0  and  R  —  R^  >  0,  subject  to  the  unconstrained  dynamics  (4). 

The  nominal  controller  is  then  augmented  with  the  nonlinear  anti-windup  scheme  described  in 
the  previous  section.  As  discussed  in  the  previous  section,  optimal  control  techniques  can  also  be 
used  to  calculate  the  gains  ki^  and  k2^.  Specifically,  the  cost  function 

is  employed  to  optimally  determine  the  gain  while  the  gain  ^2^  is  calculated  using  the  cost 
function 

j  =  t  + u^u)  dt. 
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Input 

Lower  Limit 

Upper  Limit 

5  [rad] 

0.037369 

0.3117 

ATo  [deg  R] 

34.495 

2984.5 

0.10005 

1 

Table  3.  Input  Constraints 


The  tuning  of  the  weighting  matrices  for  both  optimization  routines  is  discussed  in  the  case  study 
section.  The  choice  of  /3(x+)  for  the  air-breathing  hypersonic  vehicle  is  motivated  by  the  behavior 
of  Unfortunately,  because  we  have  multiple  inputs,  the  explicit  calculation  of  the  region  T+  is 
very  difficult.  However,  unless  the  unstable  state  reaches  a  value  which  causes  it  to  grow  without 
bound,  the  values  of  /3(x+)  can  be  set  arbitrarily  by  choosing  If  is  set  too  small  it  may  result 
in  performance  degradation  of  the  closed  loop  system.  In  the  case  of  the  air-breathing  HSV  model, 
the  unstable  state  does  not  exceed  T+  for  all  the  reference  trajectories  that  are  considered.  Thus, 
simply  setting  the  value  of  /3(x+)  to  be  large  enough  not  to  interfere  with  the  system  response  is 
sufficient  to  achieve  the  desired  response.  The  nominal  anti- windup  controller  keeps  inside 
The  function  /3(x+)  is  chosen  to  be  interpolated  from  the  data  points  /?(— 120)  =  0,  /?(— 110)  =  0, 
/?(— 100)  =  1,  /3(100)  =  1,  /3(110)  =  0,  /3(120)  =  0.  Figure  4  shows  the  Simulink  implementation  of 
the  closed  loop  system  with  anti-windup  control.  Notice  that  the  scheme  in  Figure  4  is  consistent 
with  the  one  in  Figure  2. 

IV.  Case  Study  of  Anti- Windup  Control  for  Air-breathing  HSVs 

In  this  section,  the  effectiveness  of  the  anti- windup  scheme  will  be  shown  by  means  of  simulation. 
It  will  be  demonstrated  that,  by  means  of  a  proper  tuning  of  the  weighting  matrices,  the  problems 
posed  by  input  saturation  can  be  overcome  at  least  locally.  The  tuning  of  the  weighting  matrices 
used  for  this  case  study  are  (5+  =  1  and  =  lO^^/gxg,  which  were  determined  iteratively  to  give 
favorable  results.  As  discussed  in,^  the  performance  objective  is  to  let  angle-of-attack,  velocity, 
and  altitude  to  track  step  inputs  filtered  by  three  decoupled  2^^^  order  command  shaping  filters. 
The  filtered  step  inputs  give  a  1  degree  command  change  in  angle-of-attack,  1000  ft/s  increase  in 
velocity,  and  10000  ft  change  in  altitude,  respectively.  In  this  study,  the  constraints  on  the  inputs 
are  explicitly  taken  into  account  by  the  means  of  anti- windup  control.  Though  the  input  constraints 
in  Table  3  were  chosen  perhaps  a  little  conservatively,  only  the  diffuser  area  ratio  could  really 
cause  a  problem.  The  diffuser  area  ratio  can  not  be  allowed  to  go  too  close  to  zero  because  this 
would  choke  the  flow  of  air  through  the  engine,  causing  the  engine  to  shut  down.  Thus  a  lower 
limit  of  0.1  was  chosen  to  give  the  engine  some  breathing  room. 

1.  First  Case:  Comparison  to  Nominal  Linear  Controller 

In  this  first  case,  the  effect  of  the  anti-windup  controller  is  illustrated  using  the  same  reference 
trajectories  and  nominal  controller  gains  given  in.^  The  reference  trajectories  have  settling  times 
of  3.6  [s]  for  (a,  166  [s]  for  Vr,  and  117  [s]  for  h.  This  case  provides  a  basis  of  comparison  for  the 
controller  with  and  without  anti- windup  compensation  and  input  saturation.  Figure  5(a)  shows 
the  system  response  on  the  nonlinear  model  in  the  absence  of  input  constraints.  This  is  the  same 
response  as  given  in  ^  for  the  set-point  tracking  controller.  Figure  5(c)  shows  the  nonlinear  system 
response  with  the  input  constraints  taken  into  consideration.  Notice  that  the  response  of  the 
system  with  input  constraints  is  worse  than  the  response  of  the  system  without  input  constraints. 
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Figure  4. 


Simulink  Block  Diagram  of  Nonlinear  System  with  Nominal  and  Anti-windup  Control 


particularly  for  the  output  It  is  clear  from  Figure  5(d)  the  the  inputs  do  indeed  saturate, 

which  leads  to  performance  degradation.  The  effectiveness  of  anti-windup  control  becomes  clear  in 
Figure  5(e).  The  output  response  closely  matches  the  system  response  without  input  constraints 
shown  in  Figure  5(a). 

2.  Second  Case:  Different  Nominal  Controller  Gains 

This  example  is  useful  in  illustrating  the  advantages  of  using  the  anti-windup  control  method. 
Though  the  reference  trajectories  are  still  the  same,  the  nominal  controller  gains  have  been  ad¬ 
justed  to  achieve  tighter  tracking  of  the  reference  trajectory  as  shown  in  Figure  6(a).  As  a  result, 
in  the  absence  of  input  constraints,  the  values  of  the  inputs  exceed  their  physical  limitations.  Fig¬ 
ure  6(a)  shows  the  response  of  the  linearized  system  without  constraints  on  the  inputs.  Figure  6(c) 
shows  that  with  the  input  constraints,  the  system  becomes  unstable  as  soon  as  all  three  inputs 
simultaneously  hit  their  saturations.  The  simulations  for  the  system  without  anti-windup  are  per¬ 
formed  on  the  linearized  system,  because  the  nonlinear  simulation  diverges  rather  quickly.  The 
main  advantage  of  anti-windup  control  is  clear.  In  Figure  6(c)  we  see  an  unstable  response  when 
input  constraints  are  included  in  the  system.  In  Figure  6(e)  we  see  that  the  anti-windup  control 
has  stabilized  the  system.  In  fact  the  system  response  is  very  close  to  the  response  in  the  absence 
of  input  constraints,  as  seen  in  Figure  6(a). 

3.  Third  Case:  Fast  Reference  Trajectories 

One  of  the  advantages  that  can  be  gained  through  anti-windup  control  is  the  ability  to  increase 
the  speed  of  the  reference  trajectories.  This  set  of  reference  trajectories  have  desirable  properties 
because  of  their  physically  achievable  and  desirable  settling  times  as  discussed  in.^  Simple  G-force 
calculations  show  that  a  desirable  response  should  have  a  settling  time  of  approximately  1.3  [s] 
for  cr,  102  [s]  for  Vr,  and  60  [s]  for  h.  Figures  7(a)  and  7(b)  show  the  linearized  system  response  to 
the  reference  trajectories  assuming  that  there  is  no  saturation  in  the  system,  allowing  the  inputs 
to  violate  their  constraints.  The  nonlinear  simulation  will  not  run  in  this  case  because  this  is  not 
a  physically  realizable  response.  Figures  7(c)  and  7(d)  show  the  linearized  system  response  when 
the  input  saturation  is  included,  but  the  linear  controller  is  not  modified  to  include  anti-windup 
control.  It  is  clear  from  these  graphs  that  the  response  with  input  saturation  and  without  anti¬ 
windup  control  is  unstable.  Again  this  leads  to  a  physically  impossible  situation,  so  the  nonlinear 
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simulation  will  not  run.  When  anti-windup  control  is  included  the  system  is  stabilized  even  with 
the  constraints  on  the  system.  Figure  7(e)  is  especially  encouraging,  because  it  comes  close  to 
matching  the  settling  times  of  the  reference  trajectories  with  minimal  overshoot. 

V.  Conclusions 

This  study  has  shown  that  anti- windup  redesign  is  a  viable  solution  to  the  problem  of  mitigating 
the  effect  of  input  saturation  in  tracking  control  of  exponentially  unstable  systems,  such  as  the  model 
of  air-breathing  HSV  under  investigation.  By  means  of  trial-and-error  iterations  a  good  selection 
for  the  controller  weights  can  be  found,  so  that  stable  tracking  is  achieved.  Anti-windup  redesign 
has  the  potential  to  allow  the  controller  to  achieve  faster  output  responses  with  respect  to  the 
case  in  which  the  constraints  are  implicitly  accounted  for  by  limiting  the  closed-loop  bandwidth, 
resulting  in  increased  transient  performance  and  more  robust  asymptotic  tracking. 
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Figure  5.  Case  1:  Nonlinear  System  Response 
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Figure  6.  Case  2:  Linearized  System  Response  with  and  without  Input  Constraints;  Nonlinear  System 
Response  with  Anti-windup  Control 
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Figure  7.  Case  3:  Linearized  System  Response  with  and  without  Input  Constraints;  Nonlinear  System 
Response  with  Anti-windup  Control 
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